Cotunneling mechanism of single-electron shuttling 
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The problem of electron transport by means of a dumbbell shaped shuttle in strong Coulomb 
blockade regime is solved. The electrons may be shuttled only in the cotunneling regime during 
the time spans when both shoulders of the shuttle approach the metallic banks. The conventional 
Anderson-like tunneling model is generalized for this case and the tunneling conductance is calcu- 
lated in the adiabatic regime of slow motion of the shuttle. Non-adiabatic corrections are briefly 
discussed. 
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I. INTRODUCTION 

The idea of inducing charge transport in nanode- 
vices by means of special type of time-dependent ma- 
nipulations with the parameters of a device was for- 
mulated theoretically as adiabatic pumping^ and nano- 
electromechanical shuttling^ mechanisms. Both phenom- 
ena are related to charge transfer through a nanoobject 
under perturbation V{x,t) periodic both in time t and 
a coordinate x. In case of pumping, x is not necessarily 
a spatial coordinate. It may be, e.g., the point on the 
Coulomb diamond diagram for the tunnel conductance 
G{vg^ , Vgr^) in double quantum dots, which circumscribes 
closed trajectories as a function of the gate voltages Vgi 
and Vg2 applied to two dots^"— . Not only charge but 
also spin density may be pumped within this paradigmi^ 
Electron-electron interaction plays essential part in the 
adiabatic pumping. In particular, the resonance Kondo- 
like tunneling is essential for the spin-charge separation of 
the pumping component of tunnel current (see, e.g., Ref. 
|7| for discussion of this problem and references therein to 
other theoretical investigations of adiabatic pumping in 
quantum dots). 

In pumping devices the nanoobject is immovable, and 
single-electron transport is induced by time-dependent 
control parameters. In the nanoelectromechanical shut- 
tling (NEMS) the charge and spin transport occurs due 
to an interplay of electronic and mechanical degrees of 
freedom: a nanoobject (or its parts) move periodically in 
real space between the electrodes and transport electrons 
from one electrode to another in the course of this mo- 
tion. In the "classical" NEMS^sS the shuttling instability 
develops because the charge accumulation in the shuttle 
induced by electron tunneling results in an electrome- 
chanical instability of the shuttle at a sufficiently large 
bias voltage and periodic mechanical motion of the shut- 
tle arises in the system. Mechanical motion of the shuttle 
may be also induced in the " dot" geometry, where a nano- 
island is attached to a flexible pillar or cantilever .i^^i^ As 
a rule, shuttling in these devices is substantially nonadi- 
abatic. 

Usage of flexible nanowires in the shuttling circuits 
opens new exciting possibilities. The flexural vibration 



of a suspended nanowire in combination with a scanning 
tunnel microscope working as one of the two electrodes 
may be used for realization of a charge transport in the 
shuttling regime jiiii^ On the other hand the nano-island 
may be attached to such a string. Then mechanical mo- 
tion of a shuttle will be provided by excitation of vibra- 
tional eigenmodes of this string J^ If the Coulomb block- 
ade in the island is strong enough, the single electron 
shuttling regime may be realized under certain exper- 
imental conditions. ^^ Although the shuttle motion in- 
duced by the string flexure is slow in comparison with 
characteristic tunneling rate, the adiabaticity is usually 
violated when the shuttle approaches the electrodes and 
moves away. 

In this paper we consider the single electron transport 
in suspended NEMS geometry. We study the shuttle 
shaped like a dumbbell (double quantum dot) rocking pe- 
riodically around a fixed axis. Such turnstile transports 
electrons by means of cotunneling mechanism. First, 
we formulate the problem of cyclic tunneling in general 
terms. Second, we calculate the periodic time-dependent 
conductance in the adiabatic regime. Finally, the role of 
non-adiabaticity will be discussed qualitatively. 



II. MODEL 

As is known, cotunneling mechanism is a cornerstone of 
the Kondo regime of tunnel conductance through quan- 
tum dots under strong Coulomb blockade. The cotunnel- 
ing amplitude Jir, which arises in the second order per- 
turbation theory in tunneling integrals Vi^r between the 
dot and the left and right metallic leads is exponentially 
weak, Jir ~ ViVr/A (A is the ionization energy of quan- 
tum dot). Due to formation of many-particle Abrikosov- 
Suhl resonances at the Fermi level, this amplitude ap- 
proaches unity at low temperatures (T -^ 0)r^ii^ and 
the zero bias anomaly (ZBA) in quantum conductance 
shows up in the Coulomb diamond window. ^^ 

Shuttling provides another enhancement mechanism 
for single electron tunneling, which exists even in the ab- 
sence of Kondo screening. ^^ This mechanism is an analog 
of the Debye- Waller effect in neutron scattering inten- 
sity: the average distance between the dot and the leads 



effectively reduces due to periodic motion of the shuttle 
between the source and drain electrodes. The exponen- 
tial enhancement of effective tunneling transparency is 
controlled by the ratio (x^)/xq, where the numerator is 
the average mean square displacement of the shuttle po- 
sition relative to its static distance a;o from the leads (the 
device is supposed to be symmetric). Recently we have 
described one more enhancement mechanism unrelated 
to Kondo effect, which may be realized in half-metals 
with the energy gap for minority spin carriers*^ This 
mechanism may be evinced as a finite bias anomaly in 
conductance, which arises due to opening of a resonance 
tunneling channel for the minority spin carriers. 

In this paper we consider tunneling through a dumb- 
bell shaped shuttle. A schematic sketch of the device is 
presented in Fig. [1] (left panel). We suppose that this 
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FIG. 1: Left panel: Turnstile shuttle in a form of a dumbbell 
suspended on a twisting string oriented perpendicularly to the 
tunneling plane. Right plane: Turnstile shuttle in a form of 
double dot suspended on a string vibrating in the tunneling 
plane. 



shuttle is suspended on an elastic string oriented perpen- 
dicularly to the figure plane, e.g., by means of the method 
developed in Ref. Il3|. Another option is to suspend two 
islands on a string where the standing wave polarized in 
the tunneling plane is excited (Fig. [TJ right panel). 

Let us turn to the twisting turnstile as a more simple 
in practical realization object. We consider the ideal- 
ized case where the equilibrium position with zero strain 
corresponds to the symmetric orientation of the dumb- 
bell relative to the edges. Then the torque vibration 
mode may be excited in the string, which induces the 
in-plane rotation of the dumbbell. When the rotation an- 
gle approaches 7r/2, the distances between the dots and 
both leads become small and the tunneling amplitude in- 
creases simultaneously for the source and the drain lead 
(see Fig. [2|) . If the Coulomb blockade is strong enough, 
the number of electrons in the turnstile is fixed and an 
electron from the lead can tunnel from the source into a 
turnstile only provided another electron leaves the dot si- 
multaneously. Thus effective elastic cotunneling arises as 
a ZBA in the tunnel conductance with strongly enhanced 
amplitude. 

One may say that electrons tunnel between two banks 
in a shuttling regime, where a carrier may be injected 
from a bank to the moving island during the limited time 
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FIG. 2: Cotunneling cycle of a turnstile shuttle 



span, when one of the two islands is passing this bank. In 
our idealized geometry this condition is satisfied by as- 
suming that in the equilibrium position at rotation angle 
ip = tunneling is exponentially weak and electron can 
jump onto the shuttle or leave it only at (^ — >■ ±7r/2. Sec- 
ond, the characteristic frequencies of string vibrations are 
of the order of several MHz under the realistic experimen- 
tal conditions ji^ This means that the mechanically driven 
shuttle moves slowly in comparison with the characteris- 
tic times of electron motion, and the time dependent pro- 
cesses are adiabatic for the most part of the shuttling pe- 
riod. Adiabaticity may be violated only when the shuttle 
approaches the metallic banks and tunneling amplitudes 
grow exponentially with time reaching its maximum for 
if = ±7r/2 and then again decreasing to nearly zero value. 
Thus the torque vibrations of a dumbbell-shape shuttle 
are transformed into periodically pulsing tunnel ampli- 
tudes with completely synchronized source-shuttle and 
shuttle-drain tunneling processes. Due to non-adiabatic 
character of switching, this tunneling is accompanied by 
formation of quasienergy levels, which contribute to the 
single-electron shuttling. 

To separate the single electron shuttling from Kondo 
shuttling, ^^ we consider here spinless carriers. This ap- 
proximation is valid provided the leads are fabricated 
from half-metals, i.e. materials in which the Fermi sur- 
face is formed by majority spin electrons, whereas the 
spectrum of minority spin electrons is gappediS,. The 
turnstile shuttle (TS) may be described within a frame- 



work of time-dependent Anderson model with the Hamil- 
tonian 

H = Hsh + Hi, + Hr + Htun (2.1) 

Here the first term is Hamiltonian of the moving TS 

il,h = Y. ^*''' + f^i2"i'^2; (2.2) 
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d]dj is the occupation operator for the electrons in 
the weUs of TS labeled as i = 1,2; Ei are the correspond- 
ing discrete energy levels. Usually the time dependence is 
explicitly included in the shuttle coordinate,i^^ because 
the shuttle is driven by external electric field. We suppose 
that the turnstile shuttle moves together with the acous- 
tically excited string, and the time variable appears only 
in the last term in the Hamiltonian (j2.ip (see below). The 
second term in Eq. (|2.2p is the electrostatic interaction 
between electrons located in two potential wells. There 
is no need to introduce the single site Coulomb blockade 
potential because the Pauli principle forbids double oc- 
cupation of any well. We assume that the coupling U12 
is strong enough and the TS may be occupied only by 
one electron located either in the level Ei or in the level 
E2- The electrons in the left (L) and right (R) banks are 
represented by the Hamiltonian 

Hi^+Hn^ J2 Yl '^/3««L«/^'t (2-3) 

I3=L,R K, 

with quasicontinuous spectra ej^- The last term in Eq. 
(|2.ip is the tunneling Hamiltonian 



^t- = E E [^/^.-(o«k^' + H-' 



(2.4) 



Due to periodic shuttling both wells of the TS en- 
ter tunneling contact with both banks. In the turnstile 
regime the time-dependent tunneling potential has the 
form of periodic pulses, W{t) — ^ w{t+pT2). In a rough 
approximation these pulses have a rectangular shape of 
duration Ti coming with the period T2, 



w{t) = wq [9{t + Ti/2) - 9{t - Ti/2)] . 



(2.5) 



This idealized picture of sudden approximation may be 
improved: if the real time-dependent tunnel integral for 
rotating TS is calculated, then the rectangular pulses 
transform into Gaussians, 



ui(i) = Wq exp 



t' 



[ 2(Ti/2)2 



(2.6) 



(see Fig. El . 

Electron cotunneling through moving turnstile may 
be mimicked by means of time dependent tunneling 
through an immovable two-channel double dot shown in 
Fig. 21 The appropriate time-dependence of the pa- 
rameters Wii{t) in this construction is provided by the 
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FIG. 3: (Color online) Periodic tunneling potential in sudden 
and Gaussian approximation 
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FIG. 4: Double dot with time-dependent tunneling rate. Elec- 
trostatic interaction Q prevents simultaneous occupation of 
both dots. 



time-dependent gate voltages regulating the magnitude 
of these parameters in such a way that in the first half- 
period the channels IL and 2R are open and the channels 
IR and 2i are blocked. As a result we get periodic in 
time tunnel integrals, containing two sequences of pulses 
shifted by half a period (see below) . 



III. SINGLE ELECTRON SHUTTLING IN 
PERIODIC REGIME 

We formulate the periodic time-dependent problem in 
the following way (see Fig. 2). Let the infinitesimally 
small bias voltage be directed from the left bank to the 
right bank. At some moment, say i = 0, the singly occu- 
pied shuttle with an electron captured in the site 2 is in 
the equilibrium position with the rotation angle tp = 0. 
Then within a quarter period i = T2/4 shuttle rotates to 
the angle (p — — 7r/2, the electron leaves well 2 for the 
right bank, thus releasing the Coulomb blockade C/12, so 
that the next electron from the left bank is allowed to 
enter the well 1. At t = 3T2/4 the rotation angle is 
if — 7r/2, and cotunneling process is allowed again. At 
t — T2 the shuttle loaded with one electron returns into 



initial position with (p = and an electron in the site 
2 and the next cycle begins. Two electrons are shuttled 
from the left lead to the right one during one cycle. 

This picture implies that only two charge states of 
TS are involved in cotunneling cycle, namely the one- 
electron states I A) — |10) and |A) = |01) where an elec- 
tron occupies the site 1 or 2, respectively. Two other 
states, namely the state with empty wells |A) = |00), 
and the doubly occupied state |11) appear only as vir- 
tual intermediate states. For the sake of simplicity we 
assume that the state |11) is completely suppressed by 
the Coulomb blockade potential Ui2- It is worthwhile to 
rewrite the Hamiltonian ()2.1|) in terms of the operators 
X^^ = |A)(A'|. Then the first term acquires the form 



H. 



E 

A=0,1.2 



EaX 



AA 



(3.1) 



Here notation |0) = |00), |1) = |10), |2) = |01) is used. 
We consider a symmetric shuttle, where Ei == E2 is the 
energy of a singly occupied shuttle and £^0 = is the en- 
ergy of the empty TS. The tunneling Hamiltonian reads 



Ht, 



Pk. i=l,2 



(3.2) 



Four tunneling parameters, which are assumed to be k- 
independent are grouped in two periodic series shifted by 
half period, 

WLl{t)=WR2{t)=Wa{t) 

WL2{t) = WRi{t) = Wb{t) = Wa{t- T2/2), (3.3) 
SO that 

Wait) - Y. ^it+PT2), 

P— — 00 

+ OC 

Wb{t) = Yl w'[i + (2^-1)^2/2)] (3.4) 

p— — 00 

with w{t) defined in Eq. (|2.5p or (|2.6p {p is integer). 

In order to describe the shuttling mechanism shown in 
Fig. [5] in mathematical terms, we represent the Hamilto- 
nian Htun in the form of discrete Fourier series. For this 
sake we expand the tunneling parameters p.3p : 



Wait) = J2^nC0s{n„t) (3.5) 

n=0 

00 

Wb{t) = ^(-l)"W„cos(f)„i) 



n=0 



where (J7„ = nO = 27rn/T2). The Fourier coefficients for 
two types of tunneling pulses (|2.5p . (|2.6I) are 



Wo = WqT, Wn 



2wq sin nriT 



TTO 



(3.6) 



for the rectangular pulses and 



Wo 



woT, Wn = v27ru'oTexp 



{ttut)'^ 



(3.7) 



for the Gaussian pulses. Here t = Ti/T2. 

Then inserting (|3.5p in (|3.2I) we see that the tunneling 
Hamiltonian 



TT _ rr(e) , tt{o) 



tun ' tun' 



is separated into two parts 

Hi^^ = j:^W(^Kt)aUX, + Jl.c., 



(3.8) 



(3.9) 



HZ^E.W(''\t)alX, + Ji.c. 



containing the contributions of the even and odd har- 
monics, 



W^^\t) = Y 2W2mCOs{n2mt), 



W-^Kt) = J2 2W2m+lCOs{n2ra+lt). 



(3.10) 



Even and odd creation operators for the lead and dot 
electrons are defined as 



4k = 7f ^"^i + ''««)' "°« " 71 ^"-^^ ~ "'=«^' 



(3.11) 



Xl = ^(X^o + X^O), Xl = ^(X^o - X^") 
V2 V2 



which implies 

xlx, - X 
xlx^ ^ X 



:ix' 



X' 



X12+X21), 



oo_±(^ll+^22_^12_^21)^ 



xlx^ = x"° = -(a:" - a:22 - x^^ ^ x"^^), 

XtX.^X- = i(X"-X- + X--X-)(3.12) 

The Hamiltonian of symmetric shuttle may be also 
rewritten in these variables as 



H.^^EdiX'^'^ + X'' 



(3.13) 



with Ed = El = E2. A similar transformation for the 
symmetric leads (e^L — £kR — £k) transforms the Hamil- 
tonian (|2.3p into 



ii^lcad = ^^''("eK^eK + O-l^C 



(3.14) 



As a result the Fourier transformed Hamiltonian of 
periodic shuttling is reduced to the system p.Sp . p.l3p . 
p.l4|) . Any Fourier harmonic Wn is a result of averaging 
the sequence of periodic pulses w{t) with the correspond- 
ing exponent. In the case of sudden rectangular pulses 
the convergence of Fourier harmonics p.6p is very poor, 
so that the whole series should be taken into account, 
and even in this case the problem of convergence known 
as the Gibbs defect remains. In the case of more realistic 
Gaussian pulses the Fourier harmonics p.7p fall exponen- 
tially, so the convergence of the series in the Hamiltonian 
p.8p is much better, 

Wn+2/Wn = exp[-^V2(2n + 2)] 

The parameter r may be estimated in a simplified model, 
where the torque mode is imitated by rotation of a sphere 
with the radius d along the orbit with the radius R, 
so that the distance between the leads equals 2{R + d). 
In this approximation r — d/ir^J R{d + R), so that the 
dwelling time Ti of the shuttle near the bank is controlled 
only by the geometrical factor d/R. In the torque regime 
the interval Ti is longer because the turnstile slows down 
around the turning point 



A. Time-dependent canonical transformation for 
shuttling Hamiltonian 

Having in mind the approximate adiabaticity of shut- 
tling transport, we use the method of time-dependent 
canonical transformation^^ for derivation of contunneling 
Hamiltonian as well as for calculation of current operator 
(see also Ref. l2J, where the time-independent canoni- 
cal transformation diagonalizing the Anderson impurity 
Hamiltonian in the mean-field approximation was pro- 
posed). 

As is shown in Ref. Il9l the tunneling terms may 
be ousted from the time-dependent Schrodinger opera- 
tor C = —ihd/dt + H with H = H^ + Ho defined in 
Eqs. p.8p - (I3.13|) by means of the transformation matrix 
gSc{t)+Sa{t) applied to the field operator A^ with compo- 
nents p. 111) . To eliminate the tunneling terms from the 
Shrodinger equation defined by the operator 



(P 



The coefficients Up^ are looked for in the form 



^=-'^§1 



H , 



(3.15) 



Uonit) = Uonit) + iVoi^{t). 



(3.18) 



where the functions Uenit) and Uonit) contain respec- 
tively even and odd cosine harmonics, similarly to Eq. 
p.lOp . The nonadiabatic contributions VeKit) and VoK{t) 
contain correspondingly even and odd sine harmonics. 

The procedure of canonical transformation is described 
in detail in Ref. 19. It is based on the Baker-Hausdorff 
expansion of exponential operator. The coefficients u 



(n) 

'PK 



and v^ are found from the condition of compensation 
of the term H in the expanded Hamiltonian H. Since 
the rotating turnstile is singly occupied at any moment, 
we project the effective Hamiltonian onto the subspace 
(e,o| ...|e,o). 

^ After elimination of iJtun bothj^he shuttle Hamiltonian 
Hsh and the lead Hamiltonian -fficad acquire corrections 
both time-independent and oscillating. Constant correc- 
tion to Hsh having the form of level renormalization can 
be calculated exactly, ^^ whereas the search of nonzero 
Fourier harmonics defined by Eqs. p.l5p - p.l8p is rather 
complicated procedure. 

It can be perceived that the time derivative in the 
operator £ p.lSp generates terms containing the factor 
ri/Afi;, where A^ — e^, — Ed- In what follows the terms 
VeK (t) and VoK (t) proportional to this small factor will be 
neglected. 

Then in the leading order 



Up^ = -^W^P\t) 



(3.19) 



with p = e,o. The time-dependent corrections to ii/^icad 
acquire the form of effective cotunneling Hamiltonian 



-fflcad — H 



lead 



He, 



M 



(3.20) 



These correction terms may be interpreted as a 
time-dependent generalizatio n^^i^^ of the well known 
SchriefFer- Wolff transformation!^ 



-'^cotunl^^j 



['-'KK'^eK'^eK'^ + J i^i^ia^i^ag^iX 



one should take H in the following form 



j:°,aLao«'X- + H.c.] 



(3.21) 



H^e^'+^°iH, + Ho)e 



-Sa-S, 



Si, 



Si=ih dAe^(^=+^°)(5e + 5o)e-^(^°+^') (3.16) 
Ja 



with 



^p^^ ("pkAT^V - H.c. 



(3.17) 



where 



j^^^; = Ty(p)(t)i^(^')(i)(^J- + J- 



are the components of time-dependent reflection and 
transmission amplitudes (see below). Returning back to 
the original variables for the leads L, R, we transform 



(l3:2T|) into ffcotun(t) ^ Hp + Hr with 

Hp{t)=HLL+HRR= (3.22) 



Hrit)^HLR + HRL^ (3.23) 

kk' 



When deriving this Hamiltonian, multiphcation rules 
p.l2p have been used. 



Before turning to calculation of the tunneling current, 
we simplify the effective Hamiltonian in order to retain 
only the terms responsible for the electron cotunneling 
though the turnstile. First, we omit the term Hp{t) p.22p 
containing creation of electron-hole pairs in the left and 
right leads. These terms are essential in the problems 
where the infrared orthogonality catastrophe modifies the 
current (Kondo or edge singularity mechanism). In the 
absence of such channels they yield only higher order 
corrections to the shuttling current. The remaining term 
Hrit) rewritten in the original variables reads 



+ 1 Y.iJ^" + J"')^^'^ - ^'')(«iK«««' " «L«L.') (3-24) 



r 



Since only the states near the Fermi level with £« ^ £f 
contribute to the current at small enough J7, we neglect 
the kk! dependence in the exchange integrals. 

By construction the difference J'^'^ — J°° stems from 
interlacing of the components Wa{t) and Wb{t) in the 
SW transformation. Namely, the 2m-th harmonic of this 
difference is 



jee _Joo^_^ Wa{t)Wbit) COS n2mtCOsQ2rnt, 

(3.25) 
where the bar denotes averaging over the shuttling pe- 
riod. One can check by inspection of Eqs. p.4p . (|2.6p 
and Fig. [3] that these averages are exponentially small. 
Correspondingly, the sum J^^ + J°° contains only the av- 
eraged products W^(i) = W^{t). These averages give the 
main contribution to the electron shuttling and we keep 
them in the tunneling current. Below we will be inter- 
ested only in a few first harmonics of tunneling current, 
which stem from the cotunneling parameters 



J"^ = 



8Wn^ 



m^i 



16W0W2 



cos 2J7i, 






(l-f-cos2m). 



(3.26) 



J = 1 cos fit. 



In these expressions only zero to second harmonics in Eq. 
p.lOp are taken into account. The higher orders terms 
~ cosnQt may be obtained within the same procedure. 



The three lowest harmonics of the tunneling potential 
W{t) contribute also to the diagonal part of the trans- 
formed shuttle Hamiltonian 



H. 



sh 



E'lX'^ 



E°^X°° 



(3.27) 



with 



K 

De = 4W^ + 2Wi + 8WaW2 cos 2nt 
Do = 2W^{l + cos2nt). 



(3.28) 



The periodic in time perturbation (13.281) results in 
a reconstruction of the energy spectrum of the shut- 
tle, which may be described in terms of quasienergies, 
f^ — E^ + mil (see, e.g., Ref. [2l|). These states are also 
involved in the tunneling current, but the corresponding 
tunneling channels open only in the higher than fourth 
order correction in Wpt. 



IV. ELECTRON TRANSPORT THROUGH A 
TURNSTILE SHUTTLE 

We start discussion of shuttling mechanism with gen- 
eral remarks. It should be emphasized that the TS can 
transport electrons preferably in one direction only pro- 
vided its mirror symmetries 1 -f-)- 2 and/or L ^ R are 
slightly violated. Indeed, if the state |2) with the elec- 
tron in the site 2 is chosen as the initial state, and the 
cycle shown in Fig. [5] is realized, then the electron charge 



is transported from the left to the right. If the config- 
uration |1) is chosen as the initial state, then the same 
cycle results in the electron transport from the right to 
the left. Another way to change the current direction is 
to shift to vr the phase of torque vibration. Certainly, two 
directions of single-electron cotunneling are equivalent if 
the mirror symmetry is perfect, so that the charge trans- 
port through the TS arises as a response to violation of 
this symmetry, whatever is the microscopic mechanism 
of this violation. 

To find the tunnel current through the turnstile, one 
should calculate the timc-dcpcndcnt probabilities of the 
transitions from the initial state at t = Iq to the states 
with 1,2,... 2M electron-hole pairs in the leads at i = 
MT2 in accordance with Fig. [3l In the problems of 
this type, one should distinguish between the weak cou- 
pling and strong coupling regimes. In the weak coupling 
regime the tunnel conductance is a superposition of elas- 
tic ZBA and the set of inelastic replica, in which the 
finite bias anomalies (FBA) arise because of excitations 
of higher harmonics of the periodic coupling strength (or 
higher Floquet states in terms of quasienergies) .— 1^ In 



the strong coupling regime the whole Floquet spectrum is 
involved in the time evolution of the system, determined 
by the evolution operator 



U{to,t)^\^{t)){^{t,)\. 



(4.1) 



This problem is quite complicated because the highly ex- 
cited Floquet states belonging to continuum are involved 
in tunneling and the relaxation channels, which open due 
to Auger-type processes of multiple creation of electron- 
hole pairs with finite energy and accompany the electron 
cotunneling. Leaving this regime for future studies, we 
consider here only the perturbative regime. 
The current operator is defined as 



i = e^{NR-NL) 



J 



(4.2) 



where Np — -^{J^k ^/3k^,3k)- Calculation of these deriva- 
tives by means of equations of motion gives for I following 
equation 



J 



^ - 5E(-^^^ + ^°°)(^"+^")K'«««-4.'«L. 



4- 






(4.3) 



(J- -I- J-) (Xi2 - X^') (a[,,a«, -t- a^^^,a^^) 



r 



Having in mind the exponential smallness of the terms 
oc (J*^*^ — J°°) discussed in the previous section we have 
kept in this equations only the terms corresponding to 
the processes shown in Fig. [51 where the electron trans- 
fer from the left to the right bank or vice versa is accom- 
panied by repopulation of the two "seats" in the shuttle. 
We calculate the tunneling conductance G(a;) arising 
as a response to weak periodic in time bias eV{t) by 
means of the Kubo-Greenwood formula 



Here \b) stands for the filled Fermi spheres in the leads. 
When choosing the Kubo-Greenwood equation in the 
form (|4.4p . we have taken into account the fact that the 
time-dependent terms are present both in the Hamilto- 
nian of unperturbed itinerant electrons in the two banks 
with applied bias 



Hb{t) = HR + HL-KNB. + NL)-eV{t)^^^^-^. (4.6) 



G(L.,i) 



huj 



[D{t,u)-D{t,u = 0)] (4.4) 



Here retarded Green function D{t, t') and its Fourier 
transform are defined as 



^{G\[m,i{t')]\G)Q{t-t') 
d{t' -t)D{t,t')e-"^^''-'\ (4.5) 



D{t,t') 



The initial "pure" state \G) is defined as a state with 
one electron in one of the two sites (say, 2) and the filled 
Fermi spheres in the bank electron reservoirs. This state 
is doubly degenerate: 



|G)e,o = |&)|(|l)±|2))/V2 



and in the tunneling Hamiltonian Htun- 

The leading time independent contribution to the cur- 
rent comes from the n = harmonic of the shuttling 
current (|4.3|) and is proportional given by the time- 
independent component ^ J'^^(°^/4 in the first term of 
Eq. ((O)) . By means of Eq. (|i^ and (|5^ we obtain 
the conventional expression for the zero bias conductance 



Go - ^P^l J^^(°V4|^ 



(4.7) 



where pp is the electron density of states in the leads. 
The prefactor 2 in this spinless model appeared due to the 
specific charge transfer mechanism in a turnstile shuttle: 
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unlike the standard shuttle with a single " seat" , the two- 
seat shuttle transfers two electrons within a period, as is 
seen from Fig. [2] 

Equation (|4.7p determines the time-averaged back- 
ground of a tunnel current through a turnstile. It is 
important that the coupling constant J'^'^(o) contains the 
maximal values of tunneling matrix elements wq corre- 
sponding to the closest contact between the shuttle and 
metallic banks, although the coupling strength is essen- 
tially weakened due to the factor t which characterizes 



the dwelling time of the shuttle in the nearest vicinity 
of the banks [see Eq. p.7p ]. Taking into account cor- 
rections given by the first harmonic J°"^^> we obtain the 
average conductance 



8ft 



(4.8) 



To find full time-dependent conductance let us rewrite 
Eq. (|4.3p for the current operator in the form 



m = f [/o(Xi2 + X2i)+/2(Xi^ + X^^)cos2m]^^(aLa^.-aL'«L« 



K k' 



+ |/,(Xl2 _ X^l) cos ntJ2J2 (4.'«fl« + ak'^L^ 



(4.9) 



r 



with 



In - 



h 



2W§ + wl + wi 



AWoWl + 2W1W2 



W^ + AWqW2 



(4.10) 



Substitution of p.7|) yields 



h 
h 



7r(woT)^ 

A 



[1 + 2e-y + 2e-^y\ 



^"ye-^/^[l + e-^^] 



(4.11) 



27r(u;or)2 



'\\ + 2e~y\. 



with y = (ttt)^ (see (|2.6p ). The parameter exp(— y) con- 
trols the smallness of the contribution of higher harmon- 
ics. 

Operators X^^ describing elastic processes with zero 
frequency do not contribute to the cotunneling dynam- 
ics and give the combinations X^^ ± X^^ in the current 
correlation functions in (|4.5p . After averaging these com- 
binations yield the factors 1 and 0, respectively in the 
resulting equations for the conductance. Inserting Eqs. 
(|4:9l) . ((4T0)) into Eq. (I43|) and then into Eq. g^) we 
calculate the electron-hole loops for itinerant carriers in 
the standard way and get finally 



G(0,t) 



-pI 



1 
2 



2 hh 



1 

4 



cos 2Vtt 



(4.12) 

This equation describes adiabatic two-electron shut- 
tling. We see that the oscillating motion of a turnstile 
with the frequency J7 results in an appearance of periodic 



component with the frequency 2VL in the time-dependent 
conductance. Higher harmonics 2mVL in the odd modes 
and {2m+ 1)0. in the coupling constants p.26p also may 
be taken into account. These harmonics result in an ap- 
pearance of the higher even harmonics oc cos 2nQ, with 
smaller amplitudes in oscillating conductance. An exam- 
ple of oscillating conductance including 4 harmonics is 
shown in Fig. [5] 



G(0,t) 




FIG. 5: Time dependent zero-bias anomaly in oscillating con- 
ductance G(0, t) — 5]]„ An cos 2nnt where for the sake of 
demonstration we have chosen a constant ratio An+i/An ~ 
1/2. Terms with n < 4 are taken into account. 



In principle, we may consider nonadiabatic effects as 
well. One may take into account weak non-adiabatic cor- 
rections, which stem from the terms Venit) and VoK,{t) in 
Eq. p.lSp . These corrections may be observed as phase 
shifts in the cosine functions. ^^ Truly non-adiabatic ef- 
fects arise when the contact between the shuttle and the 



banks results in reconstruction of the spectra of electrons 
in the shuttle due to a fast enough periodic passage of two 
islands forming a turnstile along the banks. Such recon- 
struction may be treated in terms of quasienergiesr^"— 
which are just a realization of the Floquet theorem for 
the time-periodic systems. In this case the energy con- 
serves in the process of cotunneling under the perturba- 
tion with the period Ti = 27r/i7 to within the " Umklapp" 
processes, Ef = Ei + mfl, where Eij are the initial and 
final states of cotunneling act. As a result, the satellites 
~ 2nhVL should arise in the conductance as finite bias 
anomalies. These satellites cannot be incorporated in 
the above calculation scheme. Development of an appro- 
priate non-perturbative description is beyond the frame- 
work of the present paper. 



V. CONCLUSIONS 

Our analysis of periodic shuttling within a time 
dependent Anderson model with periodic tunneling p.9l 



perturbing odd and even modes of the shuttle and bath 
subsystems has shown that the shuttle works as a har- 
monic analyzer, which transforms the input signal (I3.10p 
into the output adiabatic signal (|4.12l) . In the approxi- 
mate equation (I4.12p only two first harmonics Q, and 2J7 
are taken into account. Numerical estimates presented 
in Fig. [S] show that even in the case of large enough pa- 
rameter exp(— y) the higher harmonics only weakly per- 
turb the basic features of the effect, namely permanent 
background ex /q, where the magnitude of cotunneling 
strength is controlled by the parameter t'^ [see Eq. (|4.1ip , 
and the periodic temporal oscillations with the leading 
harmonic 2VL. 

We have found three realizations of this model, namely 
the turnstile suspended on an elastic string in two config- 
urations (Fig. [1} and a double quantum dot with gate- 
controlled tunneling parameters Wpi (Fig. [J]). In princi- 
ple turnstile configurations may be realized also in molec- 
ular motors. In this paper we considered only the charge 
transport induced by electrical bias, but the ratchet-like 
turnstiles could also be proposed. 
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